Coherent gluon production in very high energy heavy ion collisions 
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The early stages of a relativistic heavy-ion collision are examined in the framework of an effective 
classical SU(3) Yang-Mills theory in the transverse plane. We compute the initial energy and 
number distributions, per unit rapidity, at mid-rapidity, of gluons produced in high energy heavy 
ion collisions. We discuss the phenomenological implications of our results in light of the recent 
RHIC data. 
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The Relativistic Heavy Ion Collider (RHIC) is cur- 
rently colliding beams of gold nuclei at the highest center 
of mass energies, per nucleon, ^/snn — 200 GeV. The 
goal of these experiments is to explore strongly interact- 
ing matter, in particular the quark gluon plasma (QGP) 
predicted by lattice QCD [Q. 

The possible formation and dynamics of the QGP de- 
pends crucially on the initial conditions, namely, the 
distribution of partons in each of the nuclei before the 
collision. At high energies and for large nuclei, parton 
distributions saturate and form a color glass condensate 
(CGC). (For other pQCD based approaches, see Ref. 
and references therein.) 

The physics of saturated gluons in the CGC is as fol- 
lows. As the energies of the colliding nuclei grow (equiv- 
alently xsj *C 1), partons in the nuclear wave func- 
tions multiply until they begin to overlap in phase space. 
Repulsive interactions among the partons ensure that 
the occupation number saturates at a value proportional 
to l/a s . This phenomenon ||, is characterized by a 
bulk scale (the saturation scale) A s , where A 2 is pro- 
portional to the gluon density per unit area in a nucleon 
or nucleus. A simple saturation model for nucleons with 
A 2 = A 2 s0 (x /x Bj ) s with A s0 = 1 GeV, x = 3 • 10~ 4 
and 5 = 0.29 describes well deeply inelastic scattering 
data at the Hadron Electron Ring Accelerator (HERA) 
for XBj < 0.01 and all values Q 2 of the transverse mo- 
mentum squared from Q 2 ~ up to Q 2 = 450 GeV 2 [||. 
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In a heavy ion collision, the CGC "shatters" produc- 
ing "on shell" gluons. In this letter, we obtain non- 
perturbative expressions relating the energy and number 
distributions of produced gluons to the saturation scale 
A s of the CGC. Therefore, in principle, the saturation 
scale A s may be determined from heavy ion experiments. 

The CGC can be quantified in a classical effective field 
theory where A 2 is the only dimensionful scale [^) . When 
A 2 ^> Aq CD (for high energies and large nuclei), the 
coupling is weak: a s = a s (A 2 ) <C 1. However, the oc- 
cupation number is large, cx l/a s 3> 1. Thus weak cou- 
pling, classical methods are applicable and can be used to 



compute the classical parton distributions of nuclei |&0 . 
Recently, renormalization group methods have been de- 
veloped which systematically incorporate quantum cor- 
rections to the EFT 

The classical EFT can be applied to nuclear collisions 
|^| ^0[ . The spectrum of gluons produced when the CGC 
shatters is described by the solution of the classical Yang- 
Mills equations in the presence of two light cone sources, 
one for each nucleus, with initial conditions for the gauge 
fields given by the gauge fields of the two nuclei before 
the collision. Analytical expressions for classical gluon 
production were obtained to lowest order in the parton 
density [j3|-|l0[ . However, these are infrared divergent and 
need to be summed to all orders in the parton density. 
This was first done numerically by two of us for an SU(2) 
gauge theory jllj and non-perturbative expressions relat- 
ing the the energy Jl2| and number |ll| distributions of 
produced gluons to the saturation scale were obtained. 
Here we extend the work of Refs. 12 1^] to an SU(3) 
gauge theory 14 1. Our results can thereby be compared 
to available and forthcoming data from RHIC. 

Simulating the SU(3) theory is technically more diffi- 
cult than the SU(2) theory. For a comparable set of pa- 
rameters, the SU(3) case is about an order of magnitude 
more challenging numerically than the SU(2) one. The 
lattice formulation of the theory is described in detail in 
pT[ . The numerical techniques we use are well-known in 
lattice gauge theory, with one notable exception. Specif- 
ically, a new procedure had to be devised in order to 
determine the initial condition for the transverse compo- 
nents of the gauge fields. To this end, one must solve Eq. 
37 of Q: 



(U (1) + U (2) )(I + U t )-h.c.l =0, (1) 



where f/M^ 2 ) are SU(3) group elements corresponding 
to CGC of the two nuclei before the collision, U is the 
sought group element to the initial gauge field, and A M are 
the Gell-Mann matrices. In other words, an Hermitean 
matrix M{U) = i[(U^ + U^)(I + W) - h.c] must be 
proportional to the unit matrix. In order to solve (p]) 
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numerically for U, we first form a non-negative function 



1 



F(U) = = Tr(M ) - -(TrM) 



(2) 



If [7 satisfies (|j), attains its absolute minumum, 

F(U) = 0. Next, we minimize F(U) by relaxation. The 
relaxation equation has the form 



U = -i\^F(e n ^U)\^ =0 U 



(3) 



where 7^ are real variables, and U is the derivative with 
respect to the relaxation time t. The explicit expression 
for the right-hand side of (||) is somewhat lengthy and 
will be presented elsewhere. The relaxation equation is 
then integrated numerically to yield the initial condition 
we seek. 

In this work we will determine two observables: the 
energy and the number distribution of produced gluons. 
In doing so, we closely follow the procedure developed for 
the SU(2) case. In the continuum limit the theory con- 
tains two dimensional parameters: A s and the nuclear 
radius R. Any observable can therefore be expressed as 
a power of A s , times a function of the dimensionless prod- 
uct A S R and of the coupling constant g. 
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FIG. 1. (a) et/A s 3 as a function of tA s for A S R = 83.7. 
(b) et / A 8 A as a function of A s a for A S R = 83.7 (squares) and 
25 (circles), where a is the lattice spacing. Lines are fits of the 
form a — bx. 

For the transverse energy of gluons we get 
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^f E (A s R)Al 
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(4) 



The function Je is determined non-perturbatively as fol- 
lows. In Figure [l (a), we plot the Hamiltonian density, for 
a particular fixed value [|17| of A S R — 83.7 (on a 512 x 512 
lattice) in dimensionless units as a function of the proper 
time in dimensionless units. We note that in the SU(3) 
case, as in SU(2), er converges very rapidly to a constant 
value. The form of er is well parametrized by the func- 
tional former = a+P exp(— 77-). Here dEx/drj/TrR 2 = a 
has the proper interpretation of being the energy density 
of produced gluons, while td — l/7/A s is the "formation 
time" of the produced glue. 

In Figure [I] (b), the convergence of a to the contin- 
uum limit is shown as a function of the lattice spacing in 
dimensionless units for two values of A S R. In Ref. p2| , 
this convergence to the continuum limit was studied ex- 
tensively for very large lattices (up to 1024 x 1024 sites) 
and shown to be linear. The trend is the same for the 
SU(3) results-thus, despite being further from the con- 
tinuum limit for SU(3) (due to the significant increase 
in computer time) a linear extrapolation is justified. We 
can therefore extract the continuum value for a. We 
find / E (25) = 0.537 and / B (83.7) = 0.497. The RHIC 
value likely lies in this range of A S R. The formation 
time td = l/j/Ag is essentially the same for SU(2)-for 
A S R = 83.7, 7 = 0.362 ±0.023. As discussed in Ref. [||, 
it is - 0.3 fm for RHIC and - 0.13 fm for LHC (taking 
A s = 2 GeV and 4 GeV respectively). 

We now combine our expression in Eq. (Q) with our 
non-perturbative expression for the formation time to 
obtain a non-perturbative formula for the initial energy 
density, 
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This formula gives a rough estimate |19| ] of the initial en- 
ergy density, at a formation time of Td = ^/l I A S R where 
we have taken the average value of the slowly varying 
function 7 to be 7 = 0.34. 

To determine the gluon number per unit rapidity, we 
first compute the gluon transverse momentum distribu- 
tions. The procedure followed is identical to that de- 
scribed in Ref. ]T^] -we compute the number distribution 
in Coulomb gauge [Q, Vj_ ■ A±_ = 0. In Fig. ||(a), we plot 
the normalized gluon transverse momentum distributions 
versus &t / A s with the value A S R — 83.7, together with 
SU(2) result. Clearly, we see that the normalized re- 
sult for SU(3) is suppressed relative to the SU(2) result 
in the low momentum region. In Fig. |^(b), we plot the 
same quantity over a wider range in kx/As for two val- 
ues of A S R. At large transverse momentum, we see that 
the distributions scale exactly as N 2 — 1, the number of 
color degrees of freedom. This is as expected since at 
large transverse momentum, the modes are nearly those 
of non-interacting harmonic oscillators. At smaller mo- 
menta, the suppression is due to non-linearities, whose 
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effects, we have confirmed, are greater for larger values 
of the effective coupling A S R. 
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FIG. 2. Transverse momentum distribution of gluons, nor- 
malized to the color degrees of freedom, n (for) = /«/ (N% — I) 
(see Eq. (Q)) as a function of AsR for SU(3) (squares) and 
SU(2) (diamonds). Solid lines correspond to the fit in Eq.(Q). 

The SU(3) gluon momentum distribution can be fitted 
by the following function, 
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—f n (k T /A s ) , 



nR 2 drjdPkT g 2 
where / n (for/A s ) is 

exp Cy/£f+ m 2 /T cS 



(6) 
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a 2 A s 4 log(47rfc T /A s )fc T 4 



(k T /A s < 3) 

(k T /A s > 3) 
(7) 



with ai = 0.0295, m = 0.067A S , T cS = 0.93A S , and 
ai = 0.0343. At low momenta, the functional form is ap- 
proximately that of a Bose-Einstein distribution in two 
dimensions even though the underlying dynamics is that 
of classical fields. The functional form at high momen- 
tum is motivated by the lowest order perturbative calcu- 
lations 

Integrating our results over all momenta, we obtain for 
the gluon number per unit rapidity, the non-perturbative 
result, 
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We find that /jv(83.7) = 0.3. The results for a wide range 
of A S R vary on the order of 10% in the case of SU(2). 

The broad features of the CGC picture have recently 
been compared to the RHIC data pt^l] . We shall here 
discuss the phenomenological implications of our specific 
model in light of the recent RHIC data on multiplic- 
ity and energy distributions. The final multiplicity of 
hadrons |22l is related to the initial gluon multiplicity 
by the relation dN h /dr) = Ki ne i dN? /drj. Here Ki ne i is 
a factor accounting for 2 — > n gluon number changing 
processes which may occur at late times beyond when 
the classical approach is applicable p3fl . Moreover, if 
partial or full thermalization does occur |23|]2"4j| , the ini- 
tial transverse energy is reduced-both due to inelastic 
collisions prior to thermalization and subsequently due 
to hydrodynamic expansion- by a factor K wor k- We then 
have 
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f N (A s R) (A S R) 2 . 



(9) 



From the RHIC data at ^/snn — 130 GeV, we have 
dN h /dr)\ v = - 1000 for central collisions [§5[-f2§. For 
9 = 2 (a s = 0.33), ttR 2 = 148 fm 2 , and f N = 0.3, 
we have n ine iA 2 = 3.5 GeV 2 . Now, from Eq. (S), the 
ratio R h = dE%/dri/dN h /dri is, since Je/Jn = 1-66, 
R h = 1.66A S / K WO rk /Kind. The experimental value |2^] 
for ^/IWn = 130 GeV is R h = 0.5 GeV. Now, if we 
assume that there is no work done due to thermalization, 
K W ork = 1, we obtain from the two conditions A s = 1.02 
GeV and Ki ne i = 3.4 as the values that give agreement 
with the data. The latter value is the maximal amount of 
inelastic gluon production possible. Alternatively, if we 
assume hydrodynamic work is done, one obtains K wor k = 
( T // r j) 1 ^ 3 ; where tj and Tj are the final and initial times 
of hydrodynamic expansion respectively. This gives us 
K W ork ~ 2. Following the same analysis as previously, we 
obtain A s = 1.28 GeV and Ki ne i = 2.13. Thus, within 
the CGC approach, we are able to place bounds on both 
the saturation scale and on the amount of inelastic gluon 
production at RHIC energies. An independent method 
to extract A s directly from the data (albeit assuming 
parton-hadron duality) is to compute the relative event- 
by-event fluctuations of the gluon number |2j| . 
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